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ABSTRACT 

One ingredient in an empirical birthrate estimate for pulsar binaries is the fraction of sky subtended 
by the pulsar beam; the pulsar beaming fraction. This fraction depends on both the pulsar's opening 
angle and the misalignment angle between its spin and magnetic axes. The current estimates for 
pulsar binary birthrates are based on an average value of beaming fractions for only two pulsars, i.e. 
PSRs B1913-I-16 and B1534-I-12. In this paper we revisit the observed pulsar binaries to examine the 
sensitivity of birthrate predictions to different assumptions regarding opening angle and alignment. 
Based on empirical estimates for the relative likelihood of different beam half-opening angles and mis- 
alignment angles between the pulsar rotation and magnetic axes, we calculate an effective beaming 
correction factor, /b,eff , whose reciprocal is equivalent to the average fraction of all randomly-selected 
pulsars that point toward us. For those pulsars without any direct beam geometry constraints, we 
find that /b,eff is likely to be smaller than 6, a canonically adopted value when calculating birthrates 
of Galactic pulsar binaries. We calculate /b,off for PSRs J0737-3039A and J1141-6545, applying the 
currently available constraints for their beam geometry. As in previous estimates of the posterior 
probability density function P{TZ) for pulsar binary birthrates TZ , PSRs J0737-3039A and J1141-6545 
still significantly contribute to , if not dominate, the Galactic birthrate of tight pulsar-neutron star 
(NS) and pulsar-white dwarf (WD) binaries, respectively. Our median posterior present-day birthrate 
predictions for tight PSR-NS binaries, wide PSR-NS binaries , and tight PSR-WD binaries given a 
preferred pulsar population model and beaming geometry are 89 Myr~^, 0.5 Myr~^, and 34 Myr~^, 
respectively. For long-lived PSR-NS binaries, these estimates include a weak (xl.6) correction for 
slowly decaying star formation in the galactic disk. For pulsars with spin period between 10 ms and 
100 ms, where few measurements of misalignment and opening angle provide a sound basis for extrap- 
olation, we marginalized our posterior birthrate distribution P{TZ) over a range of plausible beaming 
correction factors. We explore several alternative beaming geometry distributions, demonstrating our 
predictions arc robust except in (untcstable) scenarios with many highly aligned recycled pulsars. 
Finally, in addition to exploring alternative beam geometries, we also briefly summarize how uncer- 
tainties in each pulsar binary's lifetime and in the pulsar luminosity distribution can be propagated 
into P{TZ). 

Subject headings: binaries: close-stars: neutron-white dwarfs-pulsars 
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1. INTRODUCTION 

Using pulsar survey selection effects to extrapolate outward to the entire Milky Way, the observed s ample of Milky 
Way f i eld binary pulsars constrains the present- d ay population and bi r thrate of these bi naries, e.g., iNaravan et al.l 
([1991 , iPhinnevI (fl99l , iCurran fc Lorimed (fT995[) , iKalogera et all (l200l , iKim et al.l (12001 , henceforth denoted KKL , 
and references t herein. Along with th e prop erties of the popul ation, this empirical b i rthrat e informs models for their 
formation, e.g., lO'Shaughnessv et al.l (|2008l ) (hereafter PSC), lO'Shaughnessv et al.l (l2009f ): detection rate estimates 
for gravitational- wave observatories like L IGO and VI RGO, e.g. Abbott et al. (2008); and even attempts to unify 
compact mergers with short 7-ray bursts (|Nakarll2007[) . Following KKL, a posterior prediction for the present-day 
birthrate (TZ) of pulsar binaries on similar evolutionary tracks to a known pulsar binary can be expressed in terms of 
the pulsar's beaming geometry (through the effective beaming correction factor /b.off), effective lifetime Toff, and the 
population distribution of individual pulsars (in luminosity and galaxy position, via Npsr)'- 

nn) = (reff/iVp,Jb,cff)7^e-(^-/^-/-.-)^ = C7^e-^^ (1) 

Summing over the individual contributions TZi from each specific pulsar binary i, a posterior prediction for the overall 
Galactic birthrate is 

viUtot) = y" n dn,v,{ni) siUtot - JI^^) • (2) 

As of 2009, the best constrained /b's for binary pulsars are available for PSRs B1913-H16 and B1534-H12 (|Kalogera et al.l 
|2001| ) . Previous wor ks taking an empirical a pproach rel i ed on these two pulsars for the beaming correction to the rate 
estimates, e.g. KKL. iKalogera "ei~al1 (l200l . 1ia^^rei~all MM \ (hereafter KKL06). The average value of /b ~ 6 based 
on PSRs B1913+16 and B1534-I-12 was used as 'canonical' value in order to calculate the birthrate (or merger rate) 
of pulsar binaries and the inferred detection rates for the gravitational-wave detectors. 

The motivation for this paper is to provide not only updated Galactic birthrates of pulsar binaries, but also to 
provide and explain more generic beaming correction factors for use in the birthrate estimates. In this work, we 
introduce an empirically-motivated beaming model, derive a probability distribution function for /b, and calculate the 
effective beaming correction factor /b,eff for two types of pulsar binaries, a pulsar with a neutron star (PSR-NS) or 
a white dwarf (PSR-WD) companion. Specific ally, we adopt cur re ntly availabl e const r aints on a misa l ignme nt angle 
a between pulsar spin and magnetic axes, e.g., IGil fc Hani ([l996h . IZha^g et all (I200I . iKoIonko et al.l (I2OOI . as well 
as the empirical relationship betwcn the half-opening angle p and pulsar spin period Ps as priors (e.g. Kramer et al. 
1998). For mildly recycled pulsars (10 ms < Pg < 100 ms), we anchor our theoretical expectations with observational 
bias; see ^for details. 

2. RATE CONSTRAINTS FOR PULSAR BINARIES 

Following KKL, we assume each pulsar i in a pulsar binary is equally visible throughout its effective lifetime 
everywhere along a fraction fb^i of all lines of sight from the pulsar, ignoring any bias introduced through time 
evolution of lu minosity and opening angleQ Consistent with observations of the Milky Way's star formation history 
(|Gilmorell200lh . we also assume the star formation rate and pulsar birthrate is constant over the past Tmw = 10 Gyr. 
Based on these two approximations, the mean number Ni of pulsars on the same evolutionary track as i that are 
visible at present can be related to its birthrate TZi by 

= f^'n,n (3) 

Conversely, Eq. ([3]) and Bayesian Poisson statistics uniquely determine our posterior prediction for each pulsar's 
birthrate shown in Eq. ([T]). 

2.1. Number of pulsars 

When inverting relation [3] to calculate a birthrate TZi, we determine Npsr following the procedure originally described 
in KKL. In order to examine the effects of newly discovered pulsars and estimated pulsar beaming fractions, we 
considered the same surveys we used in PSC. This includes all surveys listed i n KKL and three mo re surveys, i.e. the 
Swinburne inte rmediate-latitude su rvey using the Parkes multibeam system ([Edwards et al.ll200l() . the Parkes high- 
latitud e survey (jBurgav et al.ll2003l ). and the mid-latitude drift-scan survey with the Arecibo telescope lChampion et al.l 
(j2004[) . Explicitly, Npsr is the most likely number of a given pulsar binary population in our Galaxy. We estimate 
Npsr via a synthetic survey of 10^ similar pulsars pointing towards us, as described in KKL. If A^^et synthetic pulsars 
are found in our virtual survey, we estimate Npsr = ^0^/Ndeu this ratio agrees with the slope labelled a in KKL and 
defined in their Eq. (8). Within Poisson error of a few % (—1 / ^/ Ndet) , our results are consistent with a reanalysis of 
previous simulations (cf. Table 1 in PSC, from KKL06,KKL). 

2.2. Effective lifetime 

^ Pulsar spin precession occurs on a much shorter timescale and does not violate our assumption, if fi, i correctly accounts for the 
precession-enhanced extent of the emission cone. 
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In the relation 131 the effective lifetime encapulsates any and all factors needed to convert between the present-day 
number Ni and a birthrate f^^Ti- of pulsars that emit along our line of sight. For example, pulsars with a long visible 
lifetime t do not reach their equilibrium number (TZt), instead accumulating steadily with time. Assuming a steady 
birthrate, the effective lifetime for a binary contaning a pulsar i is the smaller of the Milky Way's age and the pulsars' 
lifetime: 

Ti = mill{Tmw,Tage + mm^Tmrg, Trf)) (4) 

The lifetime of a pulsar binary is defined as the sum of the current age of the system (rage) and an estimate of the 
remaining detectable lifetime. For the present age of PSRs B1913-I-16 and B1534-I-12, we use the spin-down age of a 
pulsar, an estimate b ased on extrapolating pure magnetic-dipole spindown backwards from present to some high initial 
frequencv (lArzouma nian et al. 1999a). For all other pulsars, we adopt the proper- motion corrected ages presented in 
iKiziltan fc ThorsettI (j2009l) . Because the pulsar spends most of its time at or near its current age, the error in the 
current age should be relatively small for any recycled pulsar, unless its current state is very close to the endpoint of 
recycling. For pulsar binaries considered in this work, even the most extreme possibility - a true age Tage = - changes 
lifetimes by less than 30%. The only exception is PSR J0737-3039A (see For sufficiently tight binaries (orbital 

periods less than ~ 10 hours), the remaining lifetime is usually limited by inspiral through gravitational radiation 
(|Petersl[i964f ). When the second-born pulsars are observed, however, the binary's detectable lifetime is instead limited 
by how long it radiates, as measured by the time until it reaches the pulsar "death-line" (r^) as shown in Fig. [1] 
(|Chen fc Rudermad[l99a iZhang et 311120001 [Harding et al.ll2002l: iContopoulos fc SDitkovskvll2006[ ). We note that our 
results are fairly robust to the uncertainties in the death-line: (a) merging timescales are more important in birthrate 
estimates for recycled pulsars {Tnu-g > Td), and (b) the death- line is relatively well-determined for the non-recycled 
pulsars in our sample, see, e.g., the discussion in i i3.3l of PSRs J1141-6545 and in fj3lT]of J1906-I-0746 as well as Figure 

The reconstructed birthrate for wide binaries (ti ~ 0(t,„^)) is sensitive to variations in the star formation history of 
the Milky Way. Small volumes of the Milky Way, such as the Hipparchos-scale volume of stars, can h ave 0(1) relati ve 
changes in the star-formation history fluctuation on 0(0.3—3 Gyr) timescales (see for example Fig. 4 in lGilmora ()200lD ). 
On the larger scales over which these radio pulsar surveys are sensitive, however, these fluctuations average out: see 
for ex ai nple observations of ope n clusters and well-mixed dwarf stars in Ide la Fuente Marcos fc de la Fuente Marc"oi 
(j2004[ ). iHernandez et al.l ()2001l ). and references therein. In particular, for lifetimes Tj < 0.2 Gyr relevant to the 
m ost signifi cant ti ght PSR-NS binaries, the ster form ation rate is constant to within tens of percent, from Fig. 15 
in Ide la Fuente Marcos &: de la Fuente Marco j ()2004[) . On longer timescales, observations of other disk galaxies and 
phenq n ienological models for galaxy assembly also supp ort a nearly-constant star formation rate fsee lNaab &: Ostrikeii 
(|2006[) . lSchoenrich fc BinnevI (|2009l ). iFuchs et all (|2009[ ) and references therein). For pulsar binarie s with r,; > 3 Gyr, 
observ atio ns and model s sugg est the star formation rate trends weakly upward with time; see, e.g.. lAumer &: BinnevI 
(pool and IFuchs et"all ([2009f ). 

For simplicity and to facilitate comparison with previous results, we adopt a constant star formation rate in most of 
this paper and figures; our final best estimates (Figure llip include a small correction for exponentially-decaying disk 
star formation, based on the candidate star formation history 

S*(0) ^ ' 

where t = is the present, d rawn from lAumer fc BinnevI ()200 9*): other proposals, such as profiles expected from 
a Kennicutt-Schmidt relation (|Fuchs et al.l 120091) . are easilv substituted. Assuming negligible delay between star 
formation and compact binary formation (i.e., TZ (x G), Eq. ([3]) for the average number of pulsars seen N at present in 
terms of the present-day birthrate TZ generalizes to 



N 



X. with time 



^f^'Rx{0)f Qdt = Nxn^x{Q) (6) 

J ~ mill ( Tx ■ Tm uj ) 



Based on the candidate star formation history of Eq. [51 (G) ~ 1.6 for the longest averaging time. Thus, because there 
was more star formation available to form the widest, long-lived PSR-NS binaries than if stars formed at a steady 
rate, the present-day formation rate for these wide PSR-NS binaries is roughly 1/1.6 w 0.65 times smaller than that 
shown in Fig. [8] 

2.3. Beaming Distribution 

Except for two pulsars, the empirical beaming geometry of any pulsar in a binary is not tightly constrained. Further, 
as discussed below, observations of single pulsars suggest that, even restricting to pulsars with similar evolutionary 
state (e.g., spin), the pulsar beam is randomly aligned relative to its spin axis. Nonetheless, because of the poisson 
statistics of pulsar detection (KKL), only one property of the intrinsic pulsar geometry distribution matters: the 
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We estimate the death ti mescale tj us ing Chen & Ruderman (1993) (see Table (T) [2] a nd Fig, [ll . which is typieally the shortest among 
different models presented bv lZhang et al.i (2000), Harding et aL. (,2002). .Contopoulos & Spitkovskvl ||2006|'). We note that th e uncertainty 
in for PSR J1 141-6545 is rough ly 60%, the death timescale for PSR J 1141-6545 estimated by a lChen fc RudermanI l l 19931) death-line is 
~0.10 Gyr, while IContopoulos fc Spitkovskv. 1,200^) curve predicts ~0.17 Gyr. 
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log Ps (S) 

Fig. 1. — Ps — Pb diagram superimposed with various death-lines. The stars indicate binaries whose visible lifetime is dominated by 
the death timescale: all wide PSR- NS binaries (recycled) and PSRs J1141-6545 and J190 6-I-0746 (nonrecycled) . Th e thick solid line is 
the Chen and Ruderman death-line llChen &: RudermanI 119931 ): the thin solid line is from [Contopoulos fc Spitkovskvl I I2006I) : the dotted 
red lines are from [Zhane et all 1120001) . The black solid line is the curve Pb/2Pb = lOGyr; the black dotted line is the curve where 
E = 4:n^IPsPs~^ = lO'^" erg/s, assuming a moment of inertia 1=10*^ g cm-^ IILorimer fc Kramerll2004| ). 



fraction F{Ps)[= /b off ] of all randomly selected pulsars whose beam crosses our line of sight. Assuming the pulsar 
beam drops off rapidly - typical models involve gaussian cones - the fraction of all pulsars F{Ps) of a given spin period 
and luminosity that emit towards us is well-defined and essentially independent of distance. The main ingredients 
necessary to calculate F{Ps) are the half-opening angle p of the radio beam and the misalignment angle between pulsar's 
spin and magnetic axes a. All of our analysis leading to the posterior rate prediction, Eqs. ([IJ — ([3]) generalize to an 
arbitrary beam geometry distribution upon substituting 1/ ff, ^ F. 

To calculate F, for simplicity and following historical convention we assume the bipolar pulsar beam subtends a 
hard-edged cone with opening angle p, misaligned by a < 7r/2 from the rotation axis. As the pulsar rotates, this beam 
subtends a solid angle 

/•raax(Q+p,7r/2) 

— = 27r X 2 / dcos6' , (7) 

jb Jmin(0, Q — p) 

where the beaming correction factor /b is the inverse of the fraction of solid angle the beam subtends, given open- 
ing and inclination angles (p, a). This simple model for beam geometry and its correlation with pulsar spin have 
been explore d ever since pu l sar polariza t ion d a ta and the rot a ting v e ctor model made alignment co nstraints possi- 
ble, see, e.g., iRankinI (I1993D . iGil fc Hani ()19960 . iKramer et all (|1998[) . IWeltevrede fc Johnston! (|2008[ ).In one form or 
another, this theoretically-motivated sem i-empirical correlation has b een adopted as an ing r edient in most model s 
for synthetic pulsar populations: see, e.g.. lArzoumanian et al.l (jl999aD . [Gonthier et al.l (|2004l ) [Gonthier et al.l (l2006f), 
iFaucher-Giguere fc Kaspil (|2006( ). IStorv et all (I2007D . and references therein^ Conversely, fully self-consistent popula- 
tion models which account for space distribution and kicks, luminosity evolution, beam structure and shape evolution, 
accretion during binary evolution, and spindown are required, in order to extract all the degenerate parameters that 
enter into a model by comparing its predictions with observations. 

^ Though some authors also adopt a n ori entation-dependent flux based on classical corc/cone models for pulsar emission, see 
lArzoumanian et al.l Ill999al) , IStorv et al.l 11200 71 ) and references therein, in order to illustrate the influence of spin-dependent beaming 
on rate predictions, we assume uniform emission over a cone. 
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In this paper, we focus on quantifying how much the addition of a spin-dependent opening angle p{Ps) combined 
with a misalignment angle a distribution influences birthrate estimates. As a simple approximati on, we adopt a and 
p( Ps) distributions that a re consistent with the observed pulse widths at 10% intensity level; see iZhang et alj (|200l 
and iKolonko et alJ ()2004l ). Specifically, we assume the misalignment angle is uniformly distributed between [0,2/7r]: 

V{a)da = 2/Trda . (8) 

Other distributions for a have been proposed, from a random vector V{a) ~ sin a (strongly disfavored by the high 
frequency with which low-a pulsars have been de tected; m ost detected pulsars would be orthogonal) to more tightly 
aligned distributions; see for example Fig. 3 in IKolonko e t al. (2004). While limited current observations cannot 
tightly constrain the intrinsic misalignment angle distribution, they strongly suggest tight alignments (low a) should 
be attained at least as often as a flat distribution implies. In Table [3] we compare our reference model with several 
alternative misalignment distributions; only for exceptional misalignment assumptions (e.g., p{a) cx 1/a) will our final 
results change significantly. 

For pulsars with spin period larger than 10ms, we calculate p(Pb) applying a mo del consistent with classical obser- 
vations of isolated pulsars's beam geometry (see.e.g.. iWeltevrede fc Johnstonlf2008l and references therein), as shown 
in Fig. M 

p{Ps)^5A°p-°-^ Ps > 10 ms (9a) 

The available single-pulsar data does not support a compelling model for rapidly spinning young and recycled pulsars. 
We adopt an ad-hoc power-law form as our fiducial choiceH 

= 54°(P,/10ms) P^<10ms (9b) 

In both regions, we allow for a small gaussian error in Inp to allow for model uncertainty; we conservatively adopt 
CTinp = In 1.3 (i.e., "30%" error; Figure[2]shows a 2a interval about our fiducial choice)0 Based on these two independent 
distributions, we estimate the fraction of pulsars with spin period P that point towards us as 

F{P,)^{f^')^l/^,,s (10) 

The trend of /h.n ff (fs) shown in Fig. [3] l argely agrees with previous estimates of the beaming fraction for nonrecycled 
pulsars; sec, c.g.. lTauris fc Manchester] (|1998l ). 

In addition to our fiducial choice, we have explored several other short-period beaming models. As a benchmark 
for comparison, two extreme cases are provided in Figure [5J a "very narrow" (blue) and "very wide" (green) opening 
angle model, where the P"^/^ truncates at 14° (149ms) and 54° (10ms), respectively. For relevant spin periods, these 
extreme alternatives imply noticeably different amounts of beaming correction from each other and our fiducial model, 
up to 0(x2); see Figure [H Nonetheless, our best semi-empirical estimates for binary pulsar birthrates are fairly or 
highly insensitive to the precise opening angle model p{P) adopted for rapidly spinning pulsars0 As described below, 
for pulsars with 10ms < Ps < 100ms, we anchor our theoretical bias with observed geometries of comparable binary 
pulsars. In this period interval, the choice for p{P) effectively serves as an upper bound on plausible p (and sets a 
lower bound on /b). 

Any specific choice of opening angle and misalignment distribution Via, p\Ps) determines that model's probability 
P{fb\Ps) that a randomly selected pulsar with spin period Pg has beaming correction factor 

Vifb\Ps)dfb = J 5{h - Ma, p))V{p\P,)dpV{a)da (11) 

Further, because l//ti is the probability that a given randomly selected pulsar points towards us, the distribution p,, 
of /f, among the fraction F of all pulsars aligned with our line of sight is 

^e(MPs)rfA = ^^|||p (12) 

For our fiducial beaming model. Fig. 2] compares the median value of fb for the two distributions as a function of 
p{Ps)- It shows the beaming correction factors fb for the detected population of pulsars can differ substantially from 
the underlying population, if only because those pulsars with narrow beam coverage (i.e., fb S> 1) are less likely to be 
detected. Only extremely narrow pulsar beams p < 10° will lead to a typical detected pulsar with fb ~ 6, comparable 

* The short-period extrapolation used here has no practical impact on our results; see Footnote 6. For example, we have also considered 
the much tighter beams implied by p cx P^ ; except for a handful of still-irrelevant PSR-WD binaries, our birthrate estimates are unchanged. 

Though the apparent beam size depends on the observing frequency slightly, because our standard error is typically much larger than 
the ch ange due to observing frequency, if any, we also adopt a frequency-independent emission cone; cf. the discussion in lMitra fc RankinI 
1I2OOI) and [Johnston et al. ( 20o|). 

® For example, while these alternatives can lead to slightly different values for /b.cff for tight PSR-WD binarie s, th e high birthrate and 
long period of PSR J1141-6545 makes the details of a short-period extrapolation astrophysically irrelevant; see i|3.3l The merger rate of 
tight PSR-NS binaries is dominated by PSRs J1906-I-0746 (P > 100ms) and J0737-3039A (limited both through constrained beam geometry 
and through anchored expectations from PSRs B1913-I-16 and B1534-I-12). The merger rate of wide PSR-NS binaries is dominated by two 
pulsars with P ~ 100ms, where single-pulsar observations strongly constrain reasonable p choices. 




log Ps (S) 

Fig. 2. — - Opening ang les p(Ps) versus spin period for a sample of pulsars drawn from IKramer et al.l II1998I ) and PSRs B1913+16 
HWeisberg Ik 'Iavlor1l2003) and B1534+12 HArzoumanian et al.lll993) . For pulsar binaries where multiple conflicting measurements of p 
exist, all values are shown (red points). For pulsars with Pg < 10 ms or Pg > 100 ms, we adopt an empirically-motivated relation for 
p{Ps), including p oc PT"'^ at large P; this distribution, allowing for scat ter (the shaded region) is used to calculate the effective beaming 
correction factor /b,cff over many pulsars with a similar spin; see Eq. 1101 For pulsars with spin period 10 ms < Ps < 100 ms, we define an 
alternative 'conservative' beaming model such that /(, = 6 (not shown). In this interval, to reflect our uncertainty in pulsar opening angle 
models, we assume In is uniformly distributed between the 'standard' beaming model implied by this flgure's p{P) a nd t he 'conservative' 
choice /i, = 6. Finally, the thin and dashed lines indicate two extreme alternative opening angle models discussed in H'2.'6\ 



to the measured value for knowu PSR-NS binaries that was adopted in previous analyses as the canonical value; see 
Fig. El 

2.4. Partial information and competing proposals 

In a few cases, observations Oi of pulsar i provide some information about that pulsar's geometry, such as confidence 
intervals on a and p or even a posterior distribution function V{a, p\Oi). When a unique, superior constraint exists, we 
could have simply factored this information into the average that defines F = 1 / fh,eS (Eq. 9) (e.g., a constraints for PSR 
J1141-6545; the preferred geometries of PSRs B1913+16 and 6 1534+12). However, in cases where competing prop osals 
exist (e.g., t he two models for PSR J0737-3039A suggested bv lDemorest et al.l (|2004[ ) and lFerdman et al.l (|2008[) . also 



reviewed in IKramer fc StairsI (j2008[) : see ^3.ip . or where our prediction is extremely sensitive t o unavailable priors, 
such as the mean misalignment angle of PSR Jl 141-6545, ([Manchester et al.]|2010t lKramer|[2008D ( >j3.3p . we explicitly 
provide multiple solutions in the text. When multiple competing predictions are available, our final prediction averages 
over a range of predictions between them. Explicitly, if Vf{logfh) is a distribution in X = log(/b//b,off) reflecting our 
a priori model uncertainty in /b,cff, and 7'(log72.|/b) = '^7'(7?.|/b) In 10 is our posterior estimate for the birthrate TZ 
given known beaming /b [cf. Eq. [T], then a posterior estimate that reflects uncertain beaming is 

r{\ogn)d\ogn = j dxVf{x + \ogh,,s)v{\ogn- x\f,,^,s) ■ (is) 

For the many PSR-NS binaries which have spin period between 10 ms and 100 ms,we calculate P(log R) using Eq. 
(11). Speciflcally, we average P(log R) over the beaming correction factor between /b.cff (listed in Table [T]) and the 
canonical value of 6, adopting a systematic error distribution Psysiiog fh) which is uniform in log /b between these 
limits. For pulsars with spin period outside of this range, we use Eq. (1) with /b,cff listed in Table [T] or [2] 

Similarly, if the lifetime is uncertain, one can marginalize P(log7^) over the lifetime r; if the relative likelihood of 
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log Ps (S) 

Fig. 3. — The effective beaming correction factor /b,off versus pulsar spin period Ps (black), assuming the di stribution of beaming 
geome tries described in the text. For comparison, wc show comparable estimates for the beaming correction factor from lTauris fc Manchesteil 
II1998I ) as a red curve (using their Eq. 15, but extended beyond the P = 0.1 — 5 s region they fit) as well as the obs erved Ps — /]-, combina tions 
for PSRs B1913+16 (log Pg = -1.23) and B1534+12 (log Ps = -1.42), based on the preferred val u es adopted inlKalogera et al.l 120011 ) (red 
points) and all possible combinations of values for a and p provided in lArzoumanian et al.l 1 11999 ). [Kramer et al.l 119981 ) ( red vertical bar). 
Finally, the gray shaded re gion indicates t he sp in period interval (10 ms < Ps < 100 ms), where few pulsar opening angles are available 
to constrain our model; see (Kramer et aTl ll 19981 ). Also shown in green and blue are the effective beaming factors implied for the extreme 
"narrow" and "wide" models shown in Fig. [2] 



different lifetimes Vri^ogTi) is known, then defining Y = log(T/ri), 

7'(log7^)dlog7^ = J dYVriY + logTi)V{logn + Y\n) , (14) 

Thougii usually the lifetime is relatively well determined, being dominated by relatively well determined merger or 
death timescale, we do use this expression to marginalize over the c onsiderable uncertain ty in the current age of PSR 
J0737-3039A, based on the proposed range of lifetimes presented bv lLorimer et al.l ()2007[) : see Fig. [SJ 

3. RESULTS 
3.1. Tight PSR-NS binaries 

The pulsar binaries wc consider in this work and the estimated /b,eff based on our standard model are summarized 
in Table [T] and [31 Tight PSR-NS binaries contribute to our estimate of the overall PSR-NS birthrate (which, for these 
short-lived systems, is equivalent to their merger rate). Among this set, both PSRs B1913+16 and B153 4-I-12 have 
both p and a measurements. For PSR B1913+16, we adopt fb,obs = 5.72 based on p = 12.4°(Wcisbc rg fc Tavlori 
200in and a = 156°. For PSR B1534-H2, we use fb^obs = 6.04 based on p = 4.87° and a = 114° (jArzoumanian et al.l 



19961) . [The canonical value of /b ~ 6 is obtained from the average of the fb.obs for these pulsars.] Alternative choices 
for a and p have been proposed for both pulsars (Fig. [3]) ; combining the most extreme of these options leads to values 
comparable to our model's preferred values: /b,eff = 2.2 and 1.8 for PSRs B1913-I-16 and B1534-I-12 respectively. As 
tightly beamed pulsars are unlikely to be seen in our fiducial or tightly beamed model [the solid and dashed curves in 



IKramer et al.l 119981 ) also obtained similar value for p 




p (degree) 

Fig. 4. — The median value of /t, for detectable pulsars (red; from Ve in Eq. I I12I1 ) and all pulsars (blue; from V in Eq. JTTJ) versus p, 
assuming a synthetic pulsar population with misalignment angle a uniformly distributed and beaming described by Eq. (0 ■ Also shown 
is l//b = (1 — cos p) + (7r/2 — p)sinp, an estimate of fb{p) assuming a is random; see Emmerling and Chevalier (1989). The detectable 
pulsar population will preferentially have wide beams. Only a population of pulsars with very tight opening angles arc consistent with 
the canonical value of = 6. In our "standard" model (Eqs. 7 and 8), such narrow beams are improbable for pulsars spinning with 
10ms < Ps < 100ms that is relevant to most of the known I^SR-NS binaries (see Table [ij. 



Fig. we adopt an alternative approach for pulsars with spin periods in the otherwise poorly constrained region 10 
and 100ms, the interval containing most PSR-NS binaries. 

Pulsars in tight PSR-NS binaries may have wide opening angles, given their spin periods (see Fig. [2]). If so, these 
binaries likely have /b < 6. However, the observations suggest that PSRs B1913+16 and B1534-I-12 have narrower 
beams (see Fig. [21 region defined by dotted lines). If these narrower opening angles are characteristic of tight, recycled 
pulsars in PSR-NS binaries, the appropriate /b could be closer to 6. In particular, given the similarity in spin period 
of PSR J0737-3039A to those pulsars and the lack of other constraints in that period interval (see Fig. [5]), we assume 
log/b could take on any value between log 1.5 (the value we estimate in our spin model) and log 6. 

Given posterior likelihoods, we could explicitly and systematically include observational constraints on the bea ming 
geometry of PSR J07 37-3039A as described earlier; see, e.g., the posterior constraints in lDemorest et al.1 ()2004[ ) and 
iFerdman et al.l ()2008l ). Observations support two alternate scenarios. In one, the pulse is interpreted as from a single 
highly aligned pole {a < 4°). Because of its tight alignment, in this model the beaming correction factor should 
be large: at least as large as those for binary pulsars (ft, ~ 6, assuming p ~ 30°, from p{Ps)), and potentially 
larger {ft, ~ 30 assuming p = 10°, based on observed openiiig ang les for PSRs B1913+16 and B1534-fl2). In the 
other scenario, favored by recent observations (jFerdman et al.|[2"008() . the pulse profile is interpreted as a double pole 
orthogonal rotator a ~ tt/2 with a fairly wide beam {p ~60°— 90°, consistent with p(Ps))- This latter case is consistent 
with our canonical model and leads to a comparable fb- Comparing with the assumptions presented earlier, so long as 
we ignore the possibility of tight alignment and narrow beams, our prefered model and uncertainties for PSR J0737- 
3039A already roughly incorporate its most significant modeling uncertainties. Considering that the contribution from 
PSR J1906-t-0746 is comparable with that of the PSR J0737-3039A, our best estimate for the birthrate of merging 

* As demonstrated here with p and in the conclusions with a [Table [3], the distribution of 'P(p,cc) must change dramatically to lead 
to a significant probability of detecting a pulsar with = 6. Rather than introduce strong assumptions and large systematic errors to 
enforce it, biasing our expectations about pulsars unlike PSRs B1913-|~16 and B1534+12 but similar to well-constrained isolated pulsars, 
we instead adopt a parallel approach. Our final results average between empirically-motivated theoretical priors, valid for all spin periods, 
and the assumption = 6, applied to pulsars similar to PSRs B1913+16 and B1534-I-12, with 10ms < Ps < 100ms. 
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TABLE 1 

Properties of PSR-NS binaries considered in this work. For most pulsars, fh,eff averages over the half-opening angle p and 
misalignment angle a. For PSR J0737-3039B, we adopt the preferred choice for a ~ 90° and average only over the stated uncertainties in 
p(Ps). For PSRs B19134-16 and B1534+12, where both a,p measurements are available, we adopt the values of fb,obs from lKalogera et al.l 
H2001I V The final column is C = Ti/Npsrfb, also see, Eq. {T} and Eq. (j4}. When numbers are uncertain, this table shows self-consistent 
fiducial choices. Significant uncertainties are included by explicit convolutions described in t he tex t. Small uncertainties are ignored; for 
example, our Monte Carlo estimates for Npsr have poisson sample-size errors of roughly l/v' N^et — 0{2 — 5%), where N^et = lO^/A^par- 



PSR Name 


(ms) 10-18 (3,-1) (M©) 


Mc 
(Mq) 


-^orb 

(hr) 


c 


/b.obs /b,cff 


'age 

(Gyr) 


(Gyr) 


(Gyr) 


C 
(kyr) 


Rcf'' 





tight binaries 



B1913-I-16 


59. 


8.63 


1, 


.44 


1 


.39 


7.75 


0.617 


5.72 


2.26 


0.0653 


0.301 


4.31 


576 


111 


1,2 


B1534-fl2 


37.9 


2.43 


1, 


.33 


1, 


.35 


10.1 


0.274 


6.04 


1.89 


0.200 


2.73 


9.48 


429 


1130 


3,4 


J0737-3039A 


22.7 


1.74 


1, 


.34 


1 


.25 


2.45 


0.088 




1.55 


0.142 


0.086 


14.2 


1403 


105 


5 


J0737-3039B 


2770. 


892. 










2.45 


0.088 




14. 


0.0493 




0.039 






6 


J1756-2251 


28.5 


1.02 


1, 


.4 


1 


.18 


7.67 


0.181 




1.68 


0.382 


1.65 


16.1 


664 


1821 


7 


J1906-I-0746 


144. 


20300. 


1, 


.25 


1, 


.37 


3.98 


0.085 




3.37 


0.000112 


0.308 


0.082 


192 


126 


8,9 


wide binaries 


































J1518-I-4904 


40.94 


0.028 


1, 


.56 


1, 


.05 


206.4 


0.249 




1.94 


29.2 


> TH 


51.0 


276 


18,700 


10, 11 


J1811-1736 


104.18 


0.901 


1, 


.60 


1 


.00 


451.2 


0.828 




2.92 


1.75 


> TH 


7.9 


584 


5860 


12,13 


J1829+2456 


41.01 


0.053 


1, 


.14 


1, 


.36 


28.3 


0.139 




1.94 


12.3 


> TH 


43.0 


271 


19,000 


14 


J1753-2240'= 


95.14 


0.97 


1 


.25 


1 


.25 


327.3 


0.303 




2.80 


1.4 


> TH 


8.2 


270 


13,900 


15 



^ Whenever available, we use the spin-down ages co rrected for the Shklovskii effects given in IKiziltan fc ThorscttI II2009I ). As for PSRs 
B1913+16 and B1534-I-12, we ad apt the results from IIArzoumanian et al.|[T999bf) . For PSRs J1906-I-0746, J1811-1736, J1829-I-2456, J1753- 
2240, which are not mentioned in IKiziltan fc ThorsettI I|2009I 1. we adopt the characteristic age as the current age of the pulsar. 
^ References: (1) Hulse & Talor (1975); (2) Wex, Kalogera, & Kramer (2001); (3) Wolszczan (1991); (4) Stairs et al. (2002); (5) Burgay et 
al. (2003); (6) Lync ct al. (2004); (7) Faulkner et al. (2004); (8) Lorimcr et al. (2006); (9) Kasian, and PALFA consortium (2008); (10) Nice, 
Sayer, & Taylor (1996); (11) Janssen ct al. (2008); (12) Lyne et al. (2000); (13) Kramer et al. (2003); (14) Champion et al. (2004); (15) Keith 

et al. (2008) 

The nature of the companion of PSR J1753-2240 is not yet clear, and it can be either a WD or NS l IKeith et al.ll2009l ). In this work, we 
assume PSR J1753-2240 is another wide NS-NS binary. Given that its small contribution to the total rate estimates, we note that the nature 
of the companion would not change the main results shown in this work. The masses shown for PSR J1753-2240 are half the total binary 
mass. All plausible mass pair choices lead to a merger time > lOGyr; the masses otherwise do not influence our results. 



TABLE 2 

Properties of tight PSR-WD binaries considered in this work. For all binaries here, /b,cff averages over the half-opening 
angle p and misalignment angle a. Monte Carlo sampling uncertainty in Npsr is roughly 1/ V N^^t ~ 0(3 — 5%) 



PSR Name Ps Ps Mpsr 
(ms) io-i« (=.-1) (Mq) 


Mc 


Porb 

(hr) 


e 


/b.off T-fg^ 

(Gyr) 


(Gyr) 


Td Npsr 

(Gyr) 


C Rcl^ 
(kyr) 





J075H-1807 


3.48 


0.00779 


1.26 


0.12 


6.32 


< 10-^ 


2.62 


6.66 


9.48 


> r„ 


2404 


1588 


1,2 


J1757-5322 


8.87 


0.0278 


1.35 


0.67 


10.9 


< 10-^^ 


1.26 


7.16 


8.0 


145 


1082 


7335 


3 


J1141-6545 


393.9 


4295. 


1.3 


0.986 


4.74 


0.172 


5.46 


0.00145 


0.60 


0.10 


346 


53 


4,5 


J1738+0333 


5.85 


0.0241 


1.7 


0.2 


8.5 


4 X 10~^ 


1.69 


3.71 


10.8 


> -TH 


609 


9716 


6 



For PSR J0751-I-1807, we use the spin-down ages corrected for the Shklovskii effects l IKiziltan & ThorsettI 120091 '). For other 
pulsars, we use the characteristic age. 

^ Reference: (1) Lundgren, Zepka, & Cordes (1995); (2) Nice, Stairs, & Kasian (2008) ; (3) Edwards & Bailes (2001) ; (4) Kaspi 
et al. (2000) ; (5) Bailes et al. (2003) ; (6) Jacoby (2005) 



PSR-NS binaries is not very sensitive to changes in a nearly orthogonal-rotator geometry model for PSR J0737-3039A. 
However, because we cannot rule out the most extreme scenarios for PSR J0737-3039A. for completeness wc also 
describe implications of a unipolar model: the beam shape constraints summarized by Fig. [7] translate to a prior on 
log/h that is roughly uniform between log 6 and log 30. 

In Fig. [6j we show that the probability distribution of PSR-NS merger rate with our best estimates for the beaming 
correction /b,eff, assuming the reference model of KKL06. P(7?.)'s follow from convolving together birthrate distribu- 
tions based on each individual pulsar binary, where those birthrate distributions are calculated as described in previous 
sections. Including PSR J1906+0746, we found the median PSR-NS merger rate is ~ 89 Myr~^, which is smaller than 
what we predicted in KKL06 ( ^ 123 Myr^^, cf. their peak value is 118 Myr~^), assuming the same Tage and T„irg 
listed in Tabled but used fbjmsj = 5.9, due entirely to the smaller beaming correction factors for PSR J0737-3039A 
allowed for in this work. Though our best estimate for the merger rate is slightly smaller than previous analyses, 
the difference is comparable to the Poisson-limited birthrate uncertainty and much smaller than the luminosity model 
uncertainty described in the app endix. Being nearly unchang ed, our study has astrophysical implications in agreement 
with prior work such as PSCand lO'Shaughnessv et al.l ()2009D . 

Although the merging timescale is reletively well-defined, we note that our estimate does not include 0(30%) 
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logR(yr-i) 

Fig. 5. — Solid curve: our estimate for the birthrate of pulsars like PSR J0737-3039A, based on our standard beaming geometry model 
(/b,eff = 1-55) and its spindown age (rage = 145Myr). Dotted curve: as previously, but marginalized over uncertainties in the age 
Tage- Motivated by averaging together models 1,4, and 5 in ILorimer et all 11200 71 ). assuming each model equally likely, we approximate 
age uncertainty in PSR J0737-3039A by a uniform distribution in t between 50 and 180 Myr. While significantly shorter lifetimes are 
permitted, the median birthrate for pulsars like PSR J0737-3039A increases by only 20%, from ISMyr-^o 22 Myr'^ 

uncertainties in the current binary age. In this work, for example, we fix the total age of PSR J0737-3039A to be ~230 
Myr. Fig. [5] shows a marginalized P(7?.) using the age constraints from (|Lorimer et al.llToOTi ) between 50 and 180 Myr 
(These models take into account interaction between the two pulsars; wc omit the two most extreme models 2 and 3 
with rapid magnetic field decay). 

3.2. Wide PSR-NS binaries 

In addition to those in tight orbits, some PSR-NS binaries have wide orbits (orbital period is larger than 10 hours). 
These wide binaries would never merge through gravitational radiation within a Hubble time. Their estimated death 
timescalcs imply that most of the known pulsars in wide binaries will remain visible for time comparable to or in 
excess of 10 Gyr (see Tabic [T]) and that bin aries have accumu lated over the Milky Way to their present number. In 
this work, we assume that PSR J1753-2240 (jKeith et al.ll2009[ ) is another wide PSR-NS binary. Motivated by the fact 
that spin periods of pulsars found in wide orbits are comparable to those relevant to tight PSR-NS binaries, we adopt 
the same techniques for estimating /b.off. For PSRs J1811-1736 and J1753-2240, wc adopt the /b.off values expected 
from their spin periods. On the other hand, for PSRs J1518+4904 and J1829-f 2456, we average 7'(72,|/b) over a range 
of In/b between the value predicted (In /b.off listed in Table [1]) and In 6. Our estimate for the birthrate of "wide" 
PSR-NS binaries, shown in Fig. \8\ is slightly higher than previously published estimates by PSC. The birthrate of 
wide PSR-NS binaries is not significantly changed due to the discovery of PSR J1753-2240, because of its resemblance 
to PSR J1811-17 36. Contrary to [21 the pre vious analyses had assumed that the pulsar population reached number 
equilibrium, e.g.. lO'Shaughnessv et al] (|2005[ ): our effective lifetime is 2 times shorter. At the same time, our beaming 
correction is roughly three times smaller, so the relevant factors mostly cancel; our median prediction is almost exactly 
1.5 times lower than the previous estimate. 

The reconstructed P{TV) assuming steady-state star formation has a median at 0.84 Myr~^ and is narrower than the 
previous work due to the discovery of PSR J1753-2240, with a ratio between the upper and lower rate estimate at 95% 
confidence interval, 7?.g5%/7?.5%, is 8.8 (compared to 13.8 without J1753-2240). As in PSC, the birthrate estimate of 
"tight" PSR-NS binaries is still roughly 100 times greater than that of "wide" binaries, even though both estimates have 
been modified through lifetime and /b,eff corrections. Interestingly, known "wide" PSR-NS binaries were discovered 




Fig. 6.— Distributions ^(logTe) for the birthrate of tight PSR-NS binaries: PSRs B1913+16 (red), B1534+12 (blue), J0737-3039A 
(green), J1756-2251 (orange), and J1906+0746 (black). The thick black curve indicates our best estimate for the overall Galactic birthrate 
of tight PSR-NS binaries, assuming the reference pulsar population model in KKL06; the median value is ~ 89 Myr~^. The individual 
colored curves indicate our best estimates for the birthrate of the individual binaries, based on Eq. (1) and incorporating /t, cff listed in 
Table [1] For the two PSRs with 10 ms < Ps < 100 ms and without specified beam geometry (PSRs J0737A-3039 and J175'6-2251), we 
account for uncertainty in opening-angle modeling by allowing In/b to be anywhere between our 'conservative' and 'standard' choices. 

through recycled pulsars exclusively. Following PSC, we interpret the difference in birthrate as evidence for a selection 
bias: only a few progenitors of PSR-NS binaries undergo enough mass transfer to spin up the PSR, but not enough to 
bring the orbit close enough to merge. 

3.3. Tight PSR- WD binaries 

With the excepti on of PSR J1 141-6545, the lifetimes of tight PSR- WD binaries arc limited by their gravitational 
wave inspiral time (jPetersI [19641 ). The net visible lifetimes of these binaries are significantly in excess of the age of 
the Milky Way, and arc strongly recycled with P < 10 ms. Based on comparison with measured isolated pulsars, 
a likely range for their opening angles and thus /b.cff can be estimated (see Table [5] for a summary). Combining 
this information, birthrates for the three millisecond PSR-WD binaries, PSRs J0751+1807, J1757-5322, and the 
newly discovered J1738-I-0333, are easily estimated (Fig. [TU)) . The birthrate for tight PSR-WD binaries, however, is 
dominated by the non-recycled pulsar PSR J1141-6545. For this unusual binary, the visi ble pulsar hfetime is dictated 
by its death timescale (tj), rather than its gravitational-wave merger time (r-^rg), e.g.. iKim et al.l (j2004f ). Based on 
different models for the death line, we expect the total lifetime of PSR J1141-6545 to be between between 100 — 170 
Myr; for the purposes of a rate estimate, we allow t to be logarithmically distributed within these limits. Additionally, 
recent observations imply its misalignment angle to be 160°l^g (or equivalently, 12° < a < 36°, which is adapted 
in this work) at a 68% confidence level, ([Manchester et al] 120101 : lKramerll2008() : see Fig. |9l Without this alignment 
information we would already predict fairly substantial beaming: /b,cff ~ 5.46 solely based on p{Ps) ~ 8°. 6 (Tabled)). 
With a nearly-polar beam, however, the beaming could be 10 times tighter. We therefore present two estimates for 
the birthrate of PSR J1141-6545, contingent on a: our current a priori model, shown as a solid curve in Fig. [HI and 
an estimate that includes tight beaming (dotted curve). 

Pulse profiles are an essential ingredient in empirical pulsar population modeling. T he pulse profil e of PSR J1141- 
6545 has significantly broadened, by a factor 8, since the pulsar was discovered in 2000 ()KrameHl2008D . This is mainly 
attributed to the effects of the geodetic precession. Recalculating Npgr for PSR Jl 141-6545 with the current, broader 
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Fig. 7. — Plausible beaming in PSR J0737-3039A: Solid curves are contours fh{ct, p) corresponding to = 1.1, 1.5, 2, 3.5, 6, 10, 30 (top to 
bottom). Overlaid are approximately-drawrn regions of a and p which are consistent with observations of PSR J073 7-3039A, correspondin g 
to the unipolar model (left box, small a) and two-pole model (right box); compare with contours and discussion in lDemorest et all l|2004l ). 

pulse, the most likely value of Npsr is now ^1000 instead of 350 that is estimated with the pulse profile in 2000 (see 
iKim et all ()2004D for the details) We expect that an analysis that treats both geodetic precession and the duration of 

fiulsar radio surveys imply a slightly different birthrate for the PSR J1141-6545 than wh at has be en shown in lKim et al.l 
20041 ). In this work, however, we use the pulse profile presented in the discovery paper (|Kaspi e t al. 2003) and assume 
no evolution in the pulse probile, as we only focus on the effects of the beaming correction factor with various beam 
geometry models. The effects of the pulse profile evolution of PSR J1141-6545 to the Galactic birthrate of PSR- WD 
binaries, taking into account more detailed corrections for observational biases such as the Doppler smearing, will be 
discussed in a separate paper (Kim et al. in prep). 
Our preferred birthrate estimate agrees with previously published estimates based on a "fidicual" beaming factor 

^ This preliminary result does not include up-to-date Doppler smearing, consistent with the wider pulse width. 
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Fig. 8.— Same as Fig. E] but for the wide PSR-NS population: PSRs J1811-1736 (red), J1518+4904 (blue), J1829+2456 (green), and 
J1753-2240 (orange), along with our best guess for the overall rate (thiek blaek curve). The median birthrate for v^ide PSR-NS binaries 
is 0.84 Myr~^, assuming the preferred pulsar population model in KKL06 and a steady star formation rate. The newly discovered PSR 
J1753-2240 is a cousin of PSR, J1811-1736: their contributions to the Galactic birthrate are comparable. 

/b = 6 [PSC; compare C = 53kyr here with A = 47kyr in their Table 1 for PSR J1141-6545], as well as with early 
estimates that adopt /b = 1 [Kim et aD (j2004l ) predict a peak-probability galactic birthrate of 4/f,Myr~^, based on an 
Npsr = 400] By coincidence the fiducial beaming factor agrees with our best beaming estimates based on p{P) data 
for comparable-spin pulsars. Our birthrate estimate also agrees with a recently published independent estimate of 
PSR-NS (or PSR-BH) birthrates ba sed on optical discovery of the tight WD-compact object binary SESS 1257-1-5428 
described in [Thompson et alj (|2009f ). 

3.4. Discussion 

Alternate misalignment models: In order to calculate /b,cff, we adopt a flat distribution as our standard model for 
the a distribution defined in Eq. ([8|). Several alternative distributions for pulsar misalignment have been proposed; 
some, however, do not self-consistently account for detection bias. For example, a randomly aligned beam vector on 
the sphere {p{a)da = dcosa) implies that nearl y all detected pulsar s should be nearly orthogonal rotators. The 
observation of many pulsars with smaller a, e.g.. iKolonko et ahl (j2004[ ). is not consistent with that distribution. On 
the other hand, a flat distribution and even more centrally concentrated ones (e.g., p{a) = cos ada) are plausible, 
considering observat ional biases can explain the lack of detected pulsars with small a (say, < 45° from Fig. 3 in 
IKolonko et al.l (|2004[ )). But because the expectation defining /b,cff (Eq. 9) is governed by the smallest typical beaming 
fraction of detected pulsars, the predictions for /b,cff sue fairly independent of misalignment model, assuming it is not 
too concentrated near the equator or poles. For this reason, and without more information to quantify uncertainties in 
a reconstructed a distribution with which to base a comparison (e.g., the selection bias of having an a measurement), 
we choose a fiat distribution for simplicity. 

What beaming to use?: No one beaming correction factor applies for all circumstances: even for a fixed beam ge- 
ometry distribution, different "natural" choices /b make sense for different questions [Figure H] . Beaming depends 
strongly on spin. Finally, the canonical choice /b = 6 implicitly requires exceptionally tight beaming, strong align- 

When describing preferred birthrates, we cite median probability. Earlier papers like KKL and IKim et al.l 1)20041 ) instead cite peak 
of the PDF p{R). For the lognormal and poisson distributions typical in this problem, the two disagree: the median is generally slightly 
higher. 
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Fig. 9. — Top panel: A contour plot of fb(ct, p) over the range relevant to PSR J1141-6545; contours shown are = 5,10,30 from 
top to b ottom, respectiv ely. The shade d region indicates the range of a and p implied by a preliminary analysis of recent observations 
[Manche ster et al. (20l3);|K rame 3 1 I2OOI) : inside this box /b ranges between 5 and 30. Bottom panel: 'P{TV) for PSR J1141-6545 base d on 
our standard model (/^ = 5.46; solid, also see TableO and the constraints above (averaging /b.off between 5.46 and 30 using Eq. JTSj; 
dashed). 
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Fig. 10.— Empirical distributions 'P{log7^) implied by the set of known, tight PSR-WD binaries: PSRs J1141-6545 (green), J1757-5322 
(purple), J0751+1807 (red), and J1738H-0333 (blue), given the assumptions about opening angle and misalignment. The P(7?.) for PSR 
J 1141-6545 presented here (solid green curve) is obtained by our standard prediction: assuming the pulsar is one of a family of pulsars 
with misalignment angles a uniformly distributed between and tt and p normally distributed as in Fig. [2](cf. see Fig. |9ll. Based on our 
standard beaming geometry model, PSR J1141-6545 still dominates the Galactic birthrate of PSR-WD binaries. The median rate for the 
Galactic PSR-WD binaries (thick black solid curve) is estimated to be ~ 34 Myr~^ 

ment, or both [Figure [H and Table [3] . If one number must be adopted, use the spin-dependent relation proposed by 
iTauris fc ManchesteH (|1998l ) for pulsars with P > 100ms [Figure [3]. 

4. CONCLUSIONS 

In this paper, we revise estimates for the birthrate of Galactic pulsar binaries, including new binaries as well as 
updated binary parameters and uncertainties. We describe a new quantity, the "effective beaming correction factor" 
/b.off, as a tool to permit reconstruction of Galactic binary pulsar birthrates from observations, when observations 
support not just one choice but a distributions of pulsar beam geometries. Currently, the best constrained /b's 
are available for PSRs B1913-t-16 and B1534-t-12, fb,obs ~ 6. Previous empirical birthrate estimates like KKL06 
adopted this factor for all pulsar binaries, ignoring experience from isolated pulsars' beams. Instead, we adopt random 
misalignment angles and a fiducial choice for p{Ps)- As summarized in Tables I and II, this choice produces /b,cff 
significantly smaller than the canonical value of 6 for most pulsar binaries (cf. also see Table [3] for comparisons 
between differerent beam ge ometry models), and for Pg > 10ms in agreement with the simple expression provided by 
iTauris fc ManchesteH (jl998f ). Generally, a detected pulsar should be a priori assumed to have nearly as large a beam 
(as small an /b) as its spin allows. Although possible, tight beams cross our line of sight rarely. Unless implicitly 
assumed ubiqutous, tightly beamed pulsars will not significantly impact a birthrate estimate. For pulsars with spin 
period between 10 ms and 100 ms, where few observations exist to justify extrapolations, we anchor our theoretical 
expectations for wide beams with observations of comparable pulsar binaries that suggest strong beaming, averaging 
between the two extremes. Most of the pulsar binaries that dominate PiTZ) have spin periods in intervals where p{Ps) 
is well-sampled. Of the three classes of binaries considered, only the birthrate of wide PSR-NS binaries is dominated 
by pulsars in the Icast-wcU-constrained interval Pg between 10 and 100 ms. Finally, the two binaries PSRs J1141-6545 
and J0737-3039A that dominate the birthrate of tight PSR-WD and PSR-NS binaries, respectively, have constrained 
beam geometries. In both cases, the plausible beam geometries allow a significantly narrower beam and thus a 
correspondingly higher birthrate. For these two pulsars and more generally for any pulsar binary with spin periods 
between 10 and 100 ms, empirical posterior distributions V{p,a), if made available, would help us better determine 
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TABLE 3 

Comparison of effective beaming correction factors /b,cfi predicted by different pulsar geometry models 
(columns 3 -6): F LAT (p(a) = 2/7r ; ADOPTED IN THE TEXT, SEE Eq. [SJ; p{a) = COS a; p{a) AS proposed by Eq. (1121) OF 
IZhang et al.I (120031 ) (ZJM03); and p{a) oc l/o between 3° and 90°. Reading from left to right, misalignment models 
correspond to an increasing fraction of fairly aligned PULSARS; for example, the last distribution has half of all 

BINARIES WITH a < 16°, WHICH FOR THE PULSARS OF INTEREST USUALLY CORRESPONDS TO /b > 10. NONETHELESS, THE RELEVANT 
EFFECTIVE BEAMING CORRECTION FACTOR DOES NOT INCREASE AS RAPIDLY AS THE LARGEST (OR EVEN MEDIAN) /,-,. GENERALLY, 
WHEN p(a) HAS A SUBPOPULATION OF TIGHTLY- ALIGNED PULSARS, /b,ofI IS INVERSELY PROPORTIONAL TO THE FRACTION OF 
BINARIES WITH large MISALIGNMENTS, AS WELL AS BEING DIRECTLY PROPORTIONAL TO THE MEDIAN /b FOR THE 
LARGE-MISALIGNMENT SUBPOPULATION. THEREFORE, LARGE /b,cff VALUES ARE ONLY POSSIBLE IF A SIGNIFICANT FRACTION OF 
PULSARS ARE ASSUMED TO BE TIGHTLY ALIGNED; SUCH A HYPOTHESIS IS DIFFICULT TO TEST, AS THESE TIGHTLY ALIGNED PULSARS 

WILL BE EXTREMELY DIFFICULT TO SEE. 



Name 


fb.obs 


flat 


ZJM03 


cos a 


1/a 


tight PSR-NS 












B1913-I-16 


5.72 


2.26 


2.63 


2.62 


3.37 


B1534-I-12 


6.04 


1.89 


2.14 


2.13 


2.69 


J0737-3039A 




1.55 


1.70 


1.69 


2.08 


J1756-2251 




1.68 


1.88 


1.87 


2.33 


J1906-I-0746 




3.37 


4.05 


3.98 


5.25 


wide PSR-NS 












J1518-I-4904 




1.94 


2.21 


2.21 


2.81 


J1811-1736 




2.92 


3.46 


3.42 


4.48 


J1829-I-2456 




1.94 


2.22 


2.21 


2.79 


J1753-2240 




2.8 


3.3 


3.27 


4.28 


tight PSR-WD 












J0751-I-1807 




2.62 


3.08 


3.06 


4.00 


J1757-5322 




1.26 


1.35 


1.33 


1.57 


J1141-6545 




5.46 


6.66 


6.46 


8.65 


J1738-I-0333 




1.69 


1.90 


1.89 


2.37 



model-independent detected beam geometries and therefore improve our estimates of pulsar binary birthrates. 

Studies of nonrecycled pulsars (most with P ~ 0(0.1 — 2 s)) have suggested that isol ated and binary pulsars 
could align their spin and magnetic axes (a — > 0) on 0(70, 200 Myr) timescales respectively (jWeltevrede &: JohnstonI 
l2008l) . While spin-beam alignment leads to narrower beams and therefore significantly increased typical /b for older 
nonrecycled pulsars, almost all the pulsars studied here are recycled; the two nonrec ycled PSRs J 1 906-f- 0746 and 
Jl 141-6545 are far too young for the proposed process to occur. At the other extreme, lYoung et al.I (j2009D recently 
proposed extremely rapid alignment for isolated pulsars, where on short alignment timescales ~ IMyr beams align 
and the opening angles converge to poo ~ 2 — 5°. If applied to the typical ~ O(lOOMyr) binary pulsars considered 
here, this model implies exceptionally strong beaming (/b ~ 2/p^ ~ 0(200)) and, for example, merging PSR-NS 
birthrates comparable to the galactic supernova rate. In a forthcoming paper we will address time-dependent effects 
like alignment in pulsar binaries. At present, our birthrate estimate relies only on observed pulsars as representatives 
of their evolutionary classes, without any correction fo r when along their evolutionary track they have been detected 
(e.g., in luminosity or, here, beam size; cf. lPhinnev fc'B landford (1981)). 

We find similar binary pulsar birthrates as previous studies when we adopt the same reference pulsar population 
model as KKL06 and our standard beaming geometry distribution. For example, the median birthrate of tight PSR- 
NS binaries is ~ 89 Myr^^; this estimate is lower than that of KKL06, mainly due to the narrower beams required 
there. Most pulsars in the PSR-NS population have spin period between 10 ms and 100 ms, where few measurements 
of misalignment and opening angle provide a sound basis for extrapolation. For these pulsars, we anchor our prior 
assumptions on the "fidicual" value fb~6 drawn from comparable-spin binary pulsars. We also updated the birthrate 
of wide PSR-NS binaries, including a newly discovered binary, PSR J1753-2240. Though we include improved estimates 
for the effective lifetime and beaming for each pulsar, these changes nearly cancel, leaving our birthrate estimate similar 
to PSC. The discovery of PSR J1753-2240 still improves our understanding of P(7?,), adding redundancy and reducing 
the uncertainty in the birthrate. For example, the ratio between the upper and lower limits of birthrate estimates at 
95% confidence interval is 8.8, including PSR J1753-2240: 1.5 times smaller than the previous estimate (13.8) based on 
three binaries. Fin ally, we upda ted the birthrate distribution of tight PSR-WD binaries, including the newly discovered 
PSR J1738+0333 (|Jacobvll2005h : PSR J1141-6545 stiU dominates the total birthrate. Though PSR J1141-6545's pulse 
width has evolved since its discovery, its discovery-time pulse width best characterizes how surveys found it. Until a 
future attempt at time-dependent pulsar surveys, we adopt that reference width when calculating its birthrate. 

Two pulsars which dominate their respective birthrates have weakly constrained pulsar geometries: PSR J0737- 
3039A and J1141-6545. Both pulsars could be consistent with our standard beaming model; both could admit much 
narrower beams. Information about these pulsars' geometries is not included in our preferred birthrate estimates. 
However, in the text we describe how our results change if updated beaming geometry information becomes available; 
see Figures [7] and El 

To facilitate comparison with previous results like PSC, in the text we adopted a steady star formation rate and 
did not marginalize over uncertainty in the parameters of our pulsar luminosity model. Our best estimates ncluding 
these factors arc shown in Figure 1111 see the Appendix. If future observations more tightly constrain the distribution 
p{Lmin,p) of luminosity model parameters, these final composite predictions can be easily re-evaluatcd using the 
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APPENDIX 
PULSAR LUMINOSITY FUNCTIONS 

Power law 

For clarity, in the text we discussed the impact of spin-dependent pulsar beaming on the birthrate given a fixed 
pulsar population model. The results shown in this work are based on our reference model (a power-law distribution 
for a pseudoluminosity dlogN/dlogL = —1, with 0.3 mJy kpc^ as the minimum intrinsic luminosity at 400MHz, 
Gaussian distribution in radial direction where the radial scale length is assumed to be Ro = 4kpc, and exponential 
function in z, with a scale height of Zo = 1.5kpc). See KKL, KKL06 for further details on the pulsar population 
model. Although the rate estimates are sensitive to model parameters, or the fraction of faint pulsar assumed in 
a model, the qualitative feature described here are robust with different model assumptions made. As discussed in 
KKL06, however, the reconstructed pulsar population depends sensitively on the luminosity function model. In this 
appendix, we quantify the uncertainties attributed to a pulsar luminosity function following KKL06. For example, if 
the cumulative probability of a luminosity greater than L is modeled by P(> L) = {Lmin/ LY~^ , a function with two 
parameters p and Lmm, then empirically we have found the number of pulsar binaries Np^r implied by one detection 
scales as 

iVp^^aL^^lO-i-Sf (Ala) 

Assuming that a pulsar luminosity function is similar for both millisecond pulsars {Pg < 20ms) and pulsars found in 
compact binaries, KKL06 adapt the results shown in (Cordes & C hernoff ,199*^) and estimate that the uncertainty in 
Lrnin and p can be described by the two uncorrelated distributions 

ViLrnin) = 1.22L„i„(1.7 - L„iin) irnm G [0, 1.7]mJykpc^ (Alb) 

V{p) = ^^=e-^P-^^'/^'''p o-p = 0.12 [1.4,2.6] (Ale) 



The above scaling relation and PDF allow us to generalize any result for ■p(log7?.|imi„,p) presented in the main 



text, which assumed a specific pulsar luminosity model Lmin,ref = 0.3mJykpc^ and pref ~ 2.0, to a "global" result 
Pg{logTZ) that fully marginalizes over pulsar model uncertainties: 

Pg(log7^) = j V{\og'R\ \ogL,nin,p)'P{\ogL,mn)'P{p)d\ogL„^indp (A2) 

log (5 = - \0g{Lrnin/Lmin,ref) + L6(p - Pref) (A3) 
V{\0gTZ\L„iin,p) ='P{\0gTZ - Q\L,nin,ref ,Pref) (A4) 

(A5) 

Therefore, changing variables to z = l.Q{p — Pref), we find the fully-marginalized PDF VgilogTZ) follows from the 
results presented here via convolution with an ambiguity function A: 

VgilogTZ) = V{\ogTZ\L,nin,ref ,Pref) * V {- \og L,ni„) ^^(z) (A6) 

A = r i- log L„,,n)*V{z) (A7) 

Fig. [TT] shows our estimate of the global ambiguity function, given the model of Eq. (jAl[) . This distribution is roughly 
as wide as the distributions shown in the text, whose width is limited by Poisson statistics of the number of observed 
binaries. Thus, without a better resolved luminosity model, even several new binary pulsars would not reduce the 
overall uncertainty in the rate estimate. 

Lognormal versus Power-law 

Though simple, a global power-law luminosity distribution leads to predictions th at depend sensitively on the low - 
luminosity cutoff; see Eq. lAll Recent detailed pulsar population synthesis studies bv lFaucher-Giguere &: Kaspil (j2OO60 
(henceforth FK) suggest a simpler, lognormal luminosity function fits the whole pulsar population 



p(log L)d log L = ^-<'°s L> f/2al^ LdlogL (A8) 



BINARY PULSAR BIRTHRATES WITH p{Ps) 



19 



1.4 p 
1.2 ^ 
1.0 ^ 
0.8 - 
0.6 - 

0.4 ; 

0.2 - 



0.0 '-_ 





-1.0 


-0.5 


0.0 


0.5 1.0 
X 


1.5 


2.0 


1.0 














0.8 








// \ \ 
// \ \ 

' / M 






0.6 








' / * 1 

; / \ 
/ / \ 






0.4 








' / * 
' / 
; / 


\ 

\\ 




0.2 
0.0 








/ / 
/ / 

\ / 

/ 

/ ' j/ 


\\ 
\\ 
\\ 






-8 


-7 


-6 


-5 -4 
log R (yr) 


-3 


-2 



Fig. 11. — Top panel: Rate ambiguity function A. implied by the sealing relations and luminosity model uncertainties listed in Eq. ijAl}, 
shown versus X (a dimensionless variable corresponding to logarithmic rate difFcrenccs log R/Ro). Bottom panel: our best estimates for the 
birthrate of tight PSR-NS (solid), wide PSR-NS (dotted), and tight PSR-WD binaries (dashed) including marginalization over alternative 
pulsar luminosity models. In practice, the curves shown here correspond to the thick black curves presented previously, after each has been 
smoothed with the kernel in the top panel. These results include the small star formation rate correction described in i|2.2l 
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Fig. 12. — Cumulative luminosity distributions P{> L) for (i) the canonical power-law (solid); (ii) a lognormal distribution with criog^ = 
0.9 and (logL) = —1.1 (dotted); and (iii) a lognormal distribution with (logL) = logL^iin where Lmin = 0.3mjykpc~^ is our fidicual cutoff 
for the power law distribution (dashed). 

where criog l = 0.9; see their Fig 15. As shown in Figure[T2l however, a naive lognormal distribution differs substantially 
from our canonical model: even adopting (logL) = log 0.3 (mJykpc"^), substantially larger than the best-fit model in 
FK with (logL) = —1.1, a lognormal distribution predicts more bright and fewer faint pulsars than our fiducial power- 
law luminosity model. Translating to iVpj,. and birthrates, compared to our fiducial power-law distribution roughly 
four to six times more binaries must be present in the fiducial lognormal ((logL) = —1.1) for one to be detected, 
depending on the pulsar spin and width being simulated. For comparison, for each pulsar binary our "extreme" 
lognormal distribution ((logL) = log 0.3) predicts roughly the same numbers Npsr of pulsars as the fiducial power law 
distribution^ Further, the best fit lognormal distribution from FK is really a spin-dependent luminosity L cx P'^'' P'^p , 

The lognormal distribution predicts almost all pulsars will be faint and therefore visible only from a certain characteristic distance. 
The standard power-law distribution, weighted against disk area, allows a few bright pulsars to be seen far away (i.e., P(> L)r2 oc 1). Thus 
even though the cumulative P(> L) for our ad-hoc "extreme" lognormal is always below the fiducial powerlaw, they produce comparable 
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with weakly constrained exponents ep and ep. That model implicitly introduces time-dependent selection effects, 
beyond the scope of this paper. The logormal model described above is presented as a convenient summary, not a 
fundamental distribution; FK thus don't describe how reliable its parameters are. Lacking control over the model, we 
defer a detailed discussion of lognormal luminosity models to a future paper. 
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